Boundary Wess-Zumino-Novikov-Witten Model from the Pairing Hamiltonian 
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Correlation functions in the Wess-Zumino-Novikov-Witten (WZNW) theory satisfy a system of 
Knizhnik-Zamolodchikov (KZ) equations, which involve constants of motion of an exactly solvable 
model, known as Gaudin magnet. We show that modified KZ equations, where the Gaudin operators 
are replaced by constants of motion of the closely related pairing Hamiltonian, give rise to a deformed 
WZNW model that contains terms breaking translational symmetry. This boundary WZNW model 
is identified and solved. The solution establishes a connection between the WZNW model and the 
pairing Hamiltonian in the theory of superconductivity. We also argue and demonstrate on an 
explicit example that our general approach can be used to derive exact solutions to a variety of 
dynamical systems. 

PACS numbers: 71.10.Pm, 74.20.Fg, 02.30.1k 



I. INTRODUCTION 



The Wess-Zumino-Novikov-Witten (WZNW) model 
plays an important role in physics. Historically, the 
SU(2) version of the model with topological coupling 
k = 1 was used to describe low-energy dynamics of a one- 
dimensional spin- 1/2 Heisenberg antiferromagnet. At 
higher integer couplings, k, it describes quantum crit- 
ical points in the parameter space of quantum antifcr- 
romagnetic spin-S chains, with S = k/2±. Many other 
applications of the WZNW model have emerged in vari- 
ous contexts lately (see, e.g., Ref. Q)- The Lagrangian 
formulation of the theory is given by a non-linear sigma- 
model defined in the Euclidean space by the action 3 
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where integration in the second topological Wess-Zumino 
term is over a three-dimensional ball, x = {z, z, £) G B , 
whose boundary at £ = is the two-dimensional sphere, 
S 2 = dB 3 , which corresponds to a compactified complex 
plane parametrized by (z,z) and g(z,z) € SU(2). The 
integer parameter k in Eq. ([!} is the level of the cor- 
responding conformal field theory (CFT). The WZNW 
action is invariant under conformal and non-Abelian cur- 
rent algebras. The current algebra transformations have 
a chiral structure, i.e. they act on the group element 
g(z, z) as g'(z, z) = U{z)g{z, z)U(z). Here U{z) and U{z) 
are independent elements of the group ££7(2). This prop- 
erty allows to study the holomorphic (z-dependent) and 
antiholomorphic (z-dependent) sectors of the model sep- 
arately (below, we focus on the holomorphic sector). 



The A-point correlators 
G(zi...Zn) = ((f>(zi,zi)... 



of primary fields, 

(z n ,z n )) Swznw , sat- 



isfy the Knizhnik-Zamolodchikov (KZ) equations 4 , 

G({z t })=0, (2) 
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1,1' = 1,2,... TV, and wi,i> = (zi — Zf)~ ■ Here S; = 
(£>l,S?,Sf) and Sf are SU(2) generators. Amazingly, 
operators HP in Eq. @ are formally equivalent to the 
integrals of motion of a seemingly unrelated Gaudin mag- 
net model^. The Gaudin magnet represents a quantum 
spin Hamiltonian, with effective long-range interactions 
between spins, which is exactly solvable "by design." 
Its Hamiltonian can be represented as a linear combi- 
nation of the mutually commuting integrals of motion, 

0,Wl,l' as follows: 



H G H G 



H = 2j2ziHf- 



(4) 



One general question that we formulate in this paper 
is whether it is possible to derive deformed WZNW mod- 
els, whose correlators satisfy modified KZ equations @ 
with a different set of operators Hi. Below, we answer 
this question in the affirmative by providing an exam- 
ple of this "reverse engineering approach" and finding a 
boundary WZNW model, which corresponds to the oper- 
ators Hi representing the integrals of motion of the dis- 
crete pairing Hamiltonian (Richardson model )£ ~ 13 ' 15 ~ — 
closely related to the Gaudin magnet. It descends from 
the familiar BCS Hamiltonian 
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where c] and ci s are fermion creation/annihilation op- 
erators corresponding to a single-particle state \l) with 



energies zi and spin s = ±. If A > 0, the ground state is 
a superconductor with all fcrmions paired^. Then, the 

operators c\ + c\_, Q/_C;/+, and J2s \A a cis - 1/2 J /2 be- 
come algebraically equivalent to the Pauli matrices a}, 
(T;~, and of (Anderson pseudospins) . The correspond- 
ing spin Hamiltonian is the integrable Richardson model, 
which can be presented in two identical ways 
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is the "left" boundary term, and the "right" bound- 
ary term, S^^p] is given by (j^J) with z — >• z and 
J 3 (z) — ► J 3 (z). In Eq. ©, J 3 (z) is a component of the 
"left" current in the SU(2) WZNW theory, defined in 

a standard way: J a (z) = (fc/2) tr S a g(z,z)d z g^(z,z) , 

a = ±,3. Note that due to conformal invariance the 
"left" currents do not depend on z and likewise the 
"right" currents, J a (z) do not depend on z. Note that 
the term (|9|) breaks translational invariance of the model 
and hence can be interpreted as a generalized impurity 2 . 
Below we prove that the boundary action gives rises to 
generalized KZ Eqs. ([2]) and present exact results for the 
corresponding correlation functions. 



where a 3,± are Pauli matrices, and the operators 
HP = — of/2 + AH® represent N mutually commut- 
ing Hp,Hft = conserved "currents." Note that 
J2i Hi — 0j an d hence the second term in Eq. (|7|) can 

be simplified as J2i Hf 1 — ~J2i &f/2 to give the "total 
pseudo-spin magnetization," which separates the Hilbert 
space into sectors with different numbers of Cooper pairs, 
which were actually studied in Ref. IjJ 

Note that the Gaudin model is closely related to the 
Richardson model (I6I7|) and corresponds to its infinite 
coupling limit. It is interesting to see what perturbed 
WZNW model, would correspond to the KZ equations @ 
with the operators Hi replaced with Hp. This is a key 
question addressed in this paper, but in the interest of 
practical applications, we shall consider a more general 
form of "new" operators: 



Hi [C] - - v M „„ 



U [C] d Zl U~ 1 [C] + XHf [C] , where 



U[C]HflJ- l [C}. 



(7) 



Here H[ -[C] is a rotated Gaudin Hamiltonian, with 
U[C] = cxp <J2il( z i)SfQ [C, Zi] >, C is a closed contour 

in the complex plane, q(z) is an arbitrary analytic func- 
tion inside C, and [C, z] = 1 if Z lies within the region 
enclosed by the contour and zero otherwise. We empha- 
size that Eq. ([7]) contains the conserved "currents" from 
the Richardson pairing model © in the simplest spe- 
cial case, of k = 1 (i.e., S become Pauli matrices <x/2), 
C — > Coo (i-e., the contour C^ encloses all points in C), 
and q(z) = -z/X, so that Umch = exp {- J2i ^i /(2A)}. 
We first present the main result for the boundary 

WZNW model, corresponding to operators Hi [C] defined 
in Eq. ©, 

Sbwznw [C] = Svvznw + S bound [C] + S bound [C], (8) 



where S'- 



wznw 



is the standard WZNW action (□}, 



S&undP] = ~ <**/( -)^'(-) 



(9) 



II. DERIVATION OF THE BOUNDARY WZNW 
ACTION 

We are seeking to prove that a correlation function of 
arbitrary primary fields 20 in the SU(2) boundary WZNW 
model dHJ, 

G(z 1 ,---z N ) = (<t> sl (zi)...(f> aN (z N ))s BVVZNW 

= ($ [C] 4> Sl (zi) ■ ■ ■ 4> SN (z N )) SwzNW (lO) 

satisfies the generalized KZ equations ([2]) with operators 
©. Here, $ [C] = e~ Sbmmdi - c "> and s* stands for the spin, 
0< s t < (fc/2), i = l...N. 

To solve the KZ equations @ we look for $ [C] in 
the form © $ [C] = efc dzq(z)j 3 (z) ^ where q ^ ig an an _ 
alytic and differentiable function in C, and utilize the 
two standard key ingredients of the SU(2) WZNW the- 
ory and CFTs 2 ^: (i) The crux here is the operator 
product expansion satisfied by the currents with the 
same chiralities (currents with different chiralities com- 
mute); (ii) Action of the Virasoro generators on pri- 
mary fields. Then the expression for the correlation func- 
tion G(z\, ■ ■ ■ Zn) can be simplified by contracting J z (z) 

in $[C] = J2p( l /P-)(§c dz l( z ') j3 ( z ')) P with a11 P rimar y 
fields 

G(zi, ■■■z N ) = (e^ ^S'Mc^]^ {z±) . . . ^(^)\ 



where the functional averaging with respect to the 
S'wznw is understood. Using the standard technique of 
Ref. |j and taking into account the boundary operator we 
find that indeed the following identity holds 



0= 



(/(zi)§!ep,zi]-H?[c] 



G(zi,---z N ) = 0, (11) 



where H~[C] is the rotated Gaudin Hamiltonian defined 
in ©. Note that if q(z) = -z/k\ and all Zi (i = 1 . . . N) 
are inside C, Eq. ([TT|) precisely reproduces modified KZ 
equations ([2]) with the Gaudin integrals of motion re- 
placed with those of the Richardson model with the in- 
teraction parameter A. Hence, we recover the amaz- 
ing fact that the correlation functions of the boundary 



WZNW model carry information about the exact corre- 
lation functions of the pairing model. 



III. SOLUTION OF THE 
KNIZHNIK-ZAMOLODCHIKOV EQUATIONS 

We now show that the generalized KZ Eqs. (j2j can be 
solved exactly using the standard off-shell Bethe Ansatz 

technique 21 . This method applies if all Hi are commut- 
ing, in what follows we will consider the case when all 
Zi G C. Look for the solution in an integral form 
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FIG. 1: (Color online) An example of a conformal map. Func- 
tion w — f(z) — iRi + R2 + R3 log ( |^j ) maps a disc with 
radius R to a strip of width tiR-j, centered at Im(w) — Ri. 



G(zi ■ ■ ■ z N ) 



. M 
f n du «({«a}|{Zi})V({u Q }|{z i }), (12) 
J fc=l 



where M is fixed from the condition that the correlation 
function G(z% ■ ■ ■ Zn) should be a singlet with respect to 
the global SU(2): S 3 G( Zl ■ ■ ■ z N ) = for WZNW model 4 
and also for our boundary case. The integrations are to 
be taken here over canonical cycles in the n dimensional 
complex space where points z t are excluded, with coeffi- 
cients, defined by the monodromy group of the function 
X({u a },{z t }). Then it follows that M = Yh=i s »- We 
look for eigenstates of a set of commuting Hamiltonian 
operators, Hf, of the pairing model in the form 



V(«i • • • u M \{zi}) = £t(«i) • • • SHu M )\0), (13) 



where 



N 



#(«) = £• 



si 



(14) 



and the bare vacuum state, |0), is a direct product of 
lowest weight vectors of the corresponding representa- 



tion St: SflSi 



mj|sj,mj), where rrii 



For 



example if k = 1 and all primary fields in G{z\ ■ ■ ■ zjq) 
are from s, = 1/2 representation space of SU(2), then 

M = N/2 and |0) = ( ° J ® . . . ® ( ° j . In the ba- 
sis where the primary fields are defined by spin s and 
its z-projection, m — —s . . . s, their correlation function 
(0™ 1 {z\) ■ ■ ■ 4>™£ (zn)) is connected to the general expres- 
sion (fT2")) as follows 



= (sN,mN\-"(si,mi\G(zi---zw), 



(15) 



where (si,rrii \ Sf — (si,rrii \ rrii. We found that 
x({w Q }|{z.i}) has the following form 



X ({u a }\{z t }) = Xo({u a }\{zi}) 

N 



(16) 



x exp ■ 
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A(fc + 2) 



- Xk}^ rmqjzi) +7^ u a 



where 
Xo({u a }\Ui}) = Y[(Zi-Zj)V^Y[(u a -u^Vim 
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is the known solutio n 21 i 22 to the KZ Eqs. for the canonic 
SU(2) WZNW model. In general these solutions can 
be expressed analytically in terms of multi- variable con- 
fluent hypergeometric function o 23 i 24 . Note that when 
A — >• 00, the Richardson pseudospin model reduces to the 
Gaudin magnet, and consequently \ in Eq. ()16[) reduces 
to xo- Moreover, analysis of Eqs. (fT"2|) at q(z) = —z/k\ 
suggests that the integral over u\ . . . Um hi Eq- (H2) has 
a saddle point defined by the condition 
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up 
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Interestingly, this condition coincides with Richardson 
equations for the eigenvalues of the reduced BCS Hamil- 
tonian ©. 



IV. BOSONIZED ACTION AT Jfe = 1: 
IMPLICATIONS 

The WZNW model © at k = 1 can also be realized as 
a free boson theory with central charge c = l 22 . Follow- 
ing the standard bosonization technique we introduce a 
scalar field, tp — ip(z)+<p(z), and rewrite the z-component 



of the current in the form J (z) 



/2tt 



d z cp(z). Note that 



this representation of J z is correct only locally. The full 
action ([lj at k = 1 reads 



S = —— } dzdzd z ipdz<£ H 1= 

47r7 r \J2^ 



dzzd z (p ■ 



a. c. 



,(19) 



where a. c. stands for an antiholomorphic contribution. 
We note that this bosonized version of the k — 1 action 
was discussed earlier in Refs. [lailil by Sierra. 

Here we emphasize that due to the presence of the 
boundary term defined by an arbitrary contour C and 



conformal invariance of the WZNW model make the 
boundary WZNW model a very useful tool to classify 
and study low-energy, strong coupling disordered quan- 
tum systems 25 as well as various systems driven out of 
equilibrium. Let us use bosonic version of the k = 1 the- 
ory as an illustrative example. If all Zi are real, we have 
a standard physical Richardson pseudo-spin model. Here 
contour C can be chosen as boundary of a narrow strip 
encompassing all z\, which explicitly shows that we have 
an equilibrium system. On contrary, if some of Zi have 
nonzero imaginary part, the contour can be a circle with 
radius R. In this case parametrization z = it + x, where 
t is the dimensionless time and x is the dimensionless co- 
ordinate, is inconvenient, as we generate a complicated 
time dependent term in our bosonic Hamiltonian. Inter- 
estingly enough, this, from a first sight abstract problem 
is closely related to another, well defined and physically 
motivated system. Consider the conformal map, z — > w, 
where w = it + X, which transforms a disc with radius R 
to an infinite strip, see Fig. 1. The bulk action is invari- 
ant under such transformations, while the boundary term 

will transform into - § c , dw sinh -2 ( ""^"M ip [z(w)]. 
Here contour C is the boundary of a strip of width 7ri?3 
centered at to — R\ . This term contributes to the Hamil- 
tonian and makes it time dependent. It describes a sin- 
gle instantaneous perturbation on the system at t = to, 
which however does not brake integrability. Remark- 
ably, we can extract enormous information about physi- 
cal properties of such systems by analyzing exact corre- 
lation functions. 



V. PRACTICAL APPLICATIONS OF THE 
BOUNDARY WZNW MODEL 

Boundary action Eq. © together with the expres- 
sion (|T6]) for the correlation functions represent our main 
mathematical result. As argued, it has important con- 
sequences for a variety of seemingly unrelated physi- 
cal models, notably dynamical systems. We provide 
here an explicit example of such correspondence between 
Maxwell-Bloch (MB) theory of a two-level laser—, which 
is shown to map onto the BWZNW model. Below, we 
derive for the first time an exact solution to the system 
of MB equations with damping 



d n £ + 1 £ = V 




d(P + 7 ± -p = 


N£, 


dtAT+^S-P* 


+ £" 



(20) 



-7||AA + AA , 



where £(77, £) is the complex electrical field amplitude, 
V(r),£) is the polarization of the medium, A/"(?7, £) is the 
population inversion, and rj — ilx/c and £ = Q(t — x/c) 
are given in terms of real space, x, and time, t, with Q 
being a physical constant that depends on material and 
cavity medium, and c is the speed of light. In Eq. (l2"Ul) . 
7 > is a decay rate of energy losses inside the laser 



medium and the constants 7_l and 7|| are damping coef- 
ficients of medium polarization and population inversion. 
Dissipation in the population inversion equation tends to 
return J\f to A/0/71 1 1 which is determined by the pumping. 
Amazingly, the Hamiltonian formulation of (|20[) for 
7|| = 7j_ = 7 = reduces to a set of KZ equations 
([TT]) with linear q(z) = a z 27 ' 28 . We show now that the 
BWZNW model with q(z) = az + (3z 2 and the corre- 
sponding generalized KZ Eqs. (jlip describe the system 
(I2TJ1) with finite damping parameters 7,7.1 > 0. First, 
we observe by analogy with Ref. [27j that the set of MB 
equations (I20[) with damping and pumping can be ob- 
tained from the compatibility condition of the following 
system of liner differential equations with complex spec- 
tral parameter z G C 



and 



detf> = 



Wo 



-¥) 



d z U = 



CT3 + U 



V>, 



where 



U n = - 
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-£* 



P= 2 



:^3 \ip, 



N -V 
-V* -M 



(21) 
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(23) 



According to the method of isomonodromy solutions of 
differential equations; 28 ' 29 variety of solutions of the MB 
equations can be obtained by classifying solutions of an 
auxiliary equation, d z ip(z) = A(z, £, rj)ip(z), that are con- 
sistent with the original MB equations. We found that 
consistent with (f2"Tj) choice of A(z, £, rj), which produces 
./V-soliton solutions of MB equations in the presence of 
pumping and damping, reads 



N 



M*,£,v) = tf-fcVa + X;- 



A; 



3=1 



^Nov-k 2 , 



(24) 



with parameters fc 2 , £0 £ K. Substituting expression for 
d z ip(z) together with Eq. (J2JI into Eqs. (J3TJ) , writing 
compatibility conditions and equating the residues of the 
poles at z = Zi, i = 1 . . . N, one will obtain for functions 
A*: 



dzAj = 



d v Aj — 



{zj--^ja3 + Uo,Aj 



(25) 



where 



U 



1 / N N \ 



N 



4^Vo K ~ £0)^3 - J2 



1=1 



zi 



(26) 



Eqs. (1251) admit a Hamiltonian structure with the Poisson with 
brackets 28 

{(Are)ab> {A n )cd} = 8mn((A m )ad8bc — 5 ad {A m ) bc j, (27) XMs({Ua}|{^j}) 



which corresponds to the sl(2) algebra on a chain. 
Therefore, Eqs. (f25j) acquire the form 9gAj = 
d v Aj — \Aj, H v }, with Hamiltonian op- 



{A 3 
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= E 




k=l 
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H k 


= E 

k=l 



Zk 
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2 / 



tr(A k a 3 ) 



s-s + + s + s- 



e-$> 
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,iV 



2 7 tr(A fc a 3 ) 



AT 

E 

i 



tr(A fc ^-) 

j Zfc(zfc — Zj) 



and5+'- = (E fe= i^)i2,2i. 

Quantization of the MB system implies replacement of 
Poisson brackets by commutators, {,} — > [,], and intro- 
duction of a quantum wave function, \I/(£, {zfc}). Then 
the si (2) algebra (|2"T|) acquires the following matrix real- 
ization 



A k 



si 

SI 



si 



~s k 



(29) 



which, together with transformation £ — > i£, leads to the 
set of Hamiltonian operators corresponding to (f2"U|) : 



/u 



ftfc = 
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* 2 E 
fe=i 



■2fe 



£)« + <" 



(£-£o)-2 7 z fc 



Sg-i 



€-e 



N 

2 E 



-Zfc 



(30) 



where k = 1 . . . N, and 5 = J2k S k ■ On a quantum 
level, the "wave function," ^(^, z\ . . . Zn), satisfying the 
set of "multitime," t — > (£, Zi . . . zn), Schrodinger equa- 
tions, 






h k y, 



(31) 

unambiguously determine the solution of MB equations. 
As we see, the second set of Schrodinger equations co- 
incides with generalized KZ equations (jTTJ) with the fol- 
lowing parameters: k + 2 = 1, ?(£ — £ ) = V^i 9( z ) ~ 
— z/A — 7Z 2 , and C = Coo, which encompasses all Zj — 
(2J\forj— fc 2 ) 1 / 2 , where t-,(o£K are free parameters. The 
first equation is formally an ordinary Schrodinger equa- 
tion with "time," £ = 17(£ — x/c). Therefore, this maps 
the problem onto a dynamical boundary WZNW model, 
with the boundary action Smb = § c dzq(z)J 3 (z), cor- 
relation functions of which depend on an additional pa- 
rameter, £, playing the role of time. 

By analogy with Eq. (fT2"T) . solution for <&(£, z\... zjv) 
is then found to be 

/M 
Y\_ dUkXMB({u a }\{Zi}) 
i i 



fc=l 



x V({u Q }|{*}), 



(32) 



Xo({Wa}|{^}) 
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i(£m«-M) e 2»(£-6,)E*< 



1 w — 

exp lxiE mi (^" -^0;))+E Uq J}- 



(33) 



Eqs. (|33[) and (|32[) determine the form of the wave func- 
tion of quantum states in the quantized MB system. 
Importantly, this wave function is cardinally different 
from the correlation function of primary fields in the 
bulk WZNW model. This is because the expression in 
the right-hand-side of Eq. (f3"3"f must be integrated in 
Eq. (f3"2"j) together with \0i while only the later appears 
in the WZNW model. 

In order to find the solution of the system of classi- 
cal MB equations (|3"Tj) . one should average the angular 
momentum operators, S , S 3 , with respect to the wave 
functions (|3"3"|) . Then the solution of classical MB equa- 
tions for physical quantities £, V, and J\f can be found 
from Eqs. d23), d25]), and dMl as follows: 



<]/' 



V 2 = 64A/" 2 k ( <T 



N q- N q + 

y ■ ' J y ■ 



3=1 Zj 3=1 * 



AT = 8jVo( * H 



c-&-«x; 




(34) 



These equations provide 7V-soliton solutions to the MB 
system. In general these solutions have compact integral 
representations which can be evaluated and compared 
with other numerical^ and experimental 31 data. In Ap- 
pendix we evaluate this integral for N = 2 soliton case 
and express the solution for £ , V, and J\f in terms of 
known Kummcr confluent hypergeometric functions. 



VI. CONCLUSION 

In conclusion, this work has introduced a method of 
reverse construction of boundary WZNW models from 
the generalized Knizhnik-Zamolodchikov equations satis- 
fied by the exact correlation functions and demonstrated 
the application of this method on the explicit example 
of KZ equations with conserved currents of the Richard- 
son model. Thereby, we established a direct connection 
between the discrete pairing model of superconductivity 
and the boundary WZNW model, which we identified 
and solved. We have established that the solutions of 
modified KZ equations are defined by the off-shell states 
of the Richardson model. Our construction is close in 
spirit but technically different from the BCS/CFT corre- 
spondence discussed earlier by Sierra—. Our main moti- 



vation has been to precisely identify the boundary op- 
erator in the WZNW model, which is related to the 
Richardson-type models. Our other motivation has been 
to outline a range of practical applications of the dis- 
covered correspondence, which is argued to be very wide 
and includes a variety of dynamical systems that can be 
mapped on the BWZNW theories and solved exactly in 
many cases. One such mapping and solution for the dy- 
namical system describing radiation of a two level laser 
with pumping and damping was presented. Dynamic 
properties of the laser were computed exactly exploiting 
the integrability of this latter system. 
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In order to construct the wave function (|32j) and (|33|) . we 
should integrate V(u) together with the term correspond- 
ing to xo, [(u — z\)(u — Z2)]~ 1 / 2< -- k+2 \ over u along a con- 
tour surrounding the brunch-cut at (z\,Z2) [see Eq. (|12[) 
and the discussion] . To perform the resulting integration, 
we make use of the identity 



l+a —bu 



du(u - Zl )- a {u - z 2 y 1+a e 



B(a, 1 - a)i.Fi [1 - a; 1; b(z 2 - zi)] (36) 



VII. APPENDIX 

Here we derive the analytical expression of the two- 
soliton solution of MB equations, when N = 2 and 
M = 1. For simplicity we consider the case with fi- 
nite pumping, A/"o, and medium polarization damping, 
7_i_, coefficients, but with 7 = 0. In this case we have 
one integration parameter u and two parameters, £1,22- 
Then from Eqs. {JT3J) and flUJ) it follows that 



where 1-F1 is the Kummer confluent hypergeometric func- 
tion and B(a, 1 — a) is the Beta-function^. Then for 
a = 1/2 and b = 1/A we have 



1F1 



1 1 Qi - z 2 ) 

2' ' A 



= e 2 >- Iq 



(zi - z%) 



2A 



(37) 



v(«) = E 



2 s] 



II, I) 



u — Z\ 



It, I) 



u — z 2 



\i, t). (35) 



with Iq being the modified Bessel function of zero order. 
Now, returning to the real time, £ — > — i£, by analytic 
continuation, and keeping causal behavior of the wave 
function, 'J, we will have 



$ = ^^.^j-ie-IC-fol^-ltt-foX*!-*.)! 

~i{€-£o)(zi - z 2 ) 



. 1 '(S-tQ)(»2-^l) 

(Z1-Z2) 4 e 2 I 



i(£-£o)(z2 - Zl) 



11,1) 



+ \zi -zi) *e 2 j 



II, t) 



(38) 



According to Eqs. (|34|) . to find £ one should act by the operator S^ + S 2 on the expression (|38)) for ^ and calculate 
the norm. By doing so and after some simple algebra we obtain 



£ = 



K-&I 



-|«-€o|^-|K-€o)(zi-Z2)| Ab g 



1 iU-go)U2-'l) r 

(zi - z 2 ) *e 2 7 



i(C-£o)(z 2 -zi) 



+ (22 -zi) *e 2 J 



«(£-£o)(zi -^2) 



(39) 



where the notion Abs means absolute value. Similarly, the expressions for the polarization of the medium V and the 
population inversion J\f read: 



V = 87rA/'oe- |c - &|;! ^- |K - fo)( * 1 -* a)l AbB 



I , ,_! i(S-io)(*2-zi) T 

— (Z1-Z2) 4 e 2 / 

z 2 



i(£-£o)(z2 -z x ) 



1 , ,1 <(i^iflK»ira) T 

+ — (Z2-Z1) 4 e 2 J 

zi 



*(£-£o)(zi -z 2 ) 



(40) 



and 



N = 87rAA e- | «- ?o|2 5 L - |( «-« o)(2l - Z2)l i?e 



£-& + — - — )(»- a )-*e , Jo 



*(€ — Co)(^2 — «l) 



1 1 \ i <({-t )( n -«2J 



i(£-£o)(zi -z 2 ) 



(41) 



r 



We remind the reader that in all expressions above Zj 



2A/"o?7 — &?j C? = 1)2), with constant parameters k\ 

and &2. Here, ?y = Qx/c, £ = fi(£ — x/c) with 51 being a 
physical constant that characterizes the material and the 



cavity medium. To the best of our knowledge, Eqs. (|39p . 
(|40f and (|4T|) have never been derived in the literature 
before. 
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